Note on the chemical potential of decoupled matter in the Universe 
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Textbooks on cosmology exhibit a thermodynamic inconsistency for free streaming, decoupled 
matter. It is connected here to the chemical potential, which deviates from its equilibrium value 
^ — aksT, where a is the usual parameter of the Fermi-Dirac or Bose-Einstein distribution function. 
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I. MOTIVATION 

In the textbooks on cosmology and astrophysics, such 
as the ones of Dodelson [1] , Weinberg [2] and of Giunti 
and Kim [3] , chapters on thermodynamics state that the 
pressure p and energy density u = p(? [p is the mass 
density, c the vacuum speed of light in vacuum) satisfy 
at temperature T the relation 



dp u + p 
dT ^ T ■ 



(1) 



Dodelson mentions that in principle there occurs also the 
chemical potential, adding that it is usually irrelevant in 
cosmology [4]. These authors all rush to point out that 
(1) is satisfied for photons and other relativistic parti- 
cles, where m ~ T'* and p = ^u, but they do not mention 
that it is violated in other cases. Indeed, for nonrela- 
tivistic matter, such as baryons, cold dark matter and 



non-relativisitic neutrinos, where u 



p <^ u, with 



n the particle density, this relation is violated. We seek 
here to explain the cause of this paradox. 

Our main point will be that because thermodynam- 
ics should be valid beyond equilibrium, an old thruth 
that has allowed to formulate thermodynamics for the 
glassy state [5], Eq. (10) below must substitute Eq. (1). 
Hereto we recall the relevant laws of thermodynamics, 
apply them to equilibrium and, next, to the out of equi- 
librium situation in the expanding universe. 



II. THERMODYNAMICS 

The first law states that the change of energy of the 
system, U = pc^V, is the sum of the heat added to the 
system, the work done on the system and energy gained 
by adding particles to the system. 



mechanical work for a change of the volume V , dW 
—pdV this yields 



dS = 



dU + pdV - pdN 



(3) 



The system is extensive, implying that S{XU, XV, XN) = 
XS{U, V, N). Therefore it must hold that 



5" = si{p,n)U + S2{p,n)V + S3{p,n)N, 



(4) 



for certain functions si^2,3- Considering dS from (4) and 
comparing its dU, dV and dA^ terms with (3), it is seen 
that si = 1/T, S2 — p/T, S3 — —fi/T, which are func- 
tions of p and n. These results imply the Euler relation 



U + pV-pN 

or, taken per unit volume, 

u + p — pn 
' = T • 

Consistency between (3) and (5) imposes also 



dp — ^ ^ ^ — —dT + ndp. 



(5) 



(6) 



(7) 



Inverting the relation T = T(p, n) at fixed n, one ob- 
tains p = p{T,n), which expresses all these quantities as 
functions of T and n. This allows to derive from (7) the 
relations 



dp 
df 



u + p — pn dp 

h n— — 

T dT 



(8) 



and 



dU =dQ+dW + pdN. 



(2) 



dp 
dn 



dp 
dn 



(9) 



where p is the gain in energy per particle when one keeps 
dQ =dW = 0. With the second law dQ = TdS and the 



While n is fixed for cooling a gas in a fixed, closed volume, 
we can also admit situations where n changes with T, as 



2 



it happens when coohng a gas with a movable piston or 
in the expanding universe. Then we are interested in the 
total derivative 



dp 
dT 



dp 
df 



dp 
dn 



dn 
dT' 



and a similar definition for d/i/dT. Combining (8), (9), 
(6) and (10) we obtain 

dp u + p — nu du du 
dT = ^-^+^dl = ^ + "d^' 

a result that could also have been obtained directly by 

dividing (7) by dT. This expresses Eq. (8) for the general 
case. From (5) we have the thermodynamic relation 



p — Ts 



(12) 



III. EQUILIBRIUM THERMODYNAMICS OF 
IDEAL QUANTUM GASES 

In thermal equilibrium the grand canonical distribu- 
tion function of an ideal Fermi-Dirac gas with mass m 
and chemical potential ji = aksT reads for a mode la- 
beled by q and having energy -E(q), 



^ gnc[a-/3B(q)] ^ i ga-/3iS(q) 



(13) 



n„=0 



where /? = l/ksT. Accounting for all modes yields 



Z = [] (l + e"-''^(i)) 



(14) 



The Fermi-Dirac distribution reads 



/(q) 



1 



(15) 



We are interested in an ideal gas of particle with mass 
m and momentum q. Taking periodic boundary condi- 
tions for a cube with size V^^^, the allowed momenta are 
{QxjQyyQz) = 27r/iF~^/'^(na;, rij^, n^), with integer values 
of nx,y,z, so that 



E{q) = ^m^d^^q^c^, 



(16) 



and 



log 



Prom this we derive 



d log S = ^ log 2: Nda - [/d/3, (18) 



(10) in which we may identify U = Vu and N = Vn from 



d^q 
{2Trhy 



■f, 



I 



d^q 



and conclude that 



p = nksT log Z. 



(19) 



(20) 



(21) 



Reinserting this into (18) we verify the central relation 
(11). We finally find from (6) that the entropy density 
equals 



s = ks 



I 



^[-/log/-(l-/)log(l-/)], (22) 



To derive this result we may use that 



(3u — an = J 
and 

pp 



d^qf{pE-a) 



{2'Khf 



r _d^ 
J W 



d\ 



■/log 



1-/ 



^''^ log(l + e«-^^) 



/ 



(27r/i)3 
d^q 



log- 



(27r/l)3 ^1-/- 



(23) 



For later use wc point out that by inserting 
J2i^idqi/dqi = 3 in the first identity and performing 
a partial integration, we obtain the equivalent results 



d^q J. q^c^ 



{2-Khf ' 3E ■ 



(24) 



We can also verify from (12) that a takes the equilib- 
rium value 

a = Pfj,. (25) 
For bosons one has the grand canonical partition sum 



(26) 



the Bose-Einstein distribution 



3 



/(q) = 



1 



g/3£;(q)-a _ I ■ 



(27) 



while the entropy reads 



s — kf 



■[-/log/ + (l + /)log(l + /)]. (28) 



IV. DECOUPLED QUANTUM MATTER IN 
THE EXPANDING UNIVERSE 



With / from (34) or, more generally from (31), the 
entropy density (28) is continuous at Ti and scales as 
after that, as does the particle density (19), so that the 
entropy per particle 



(7 = 



n{t) n{Ti) 



(35) 



is constant, and the entropy in a comoving unit volume 
S = Na ^ a^s conserved. The same holds for the co- 
moving particle number N ~ a'^n. 
The energy density is thus simply 



In the expanding universe a certain species may decou- 
ple from the other matter, meaning that, given its cross 
section, the scattering candidates become so sparse that 
practically no scattering will take place anymore. Then 
in a flat Friedman metric 



ds^ = c^de - a2(t)(da;2 + dt/^ + d^^) 



(29) 



the one-particle occupation of the mode with momentum 
q satisfies the free-streaming Boltzmann equation [1-3] 



dt}{q,t)-\-^}{q,t) 
a aq 



0, 



(30) 



with right hand side equal to zero (absence of scattering), 
and a the scale factor of (29). The solution reads 



f{q,t) = f[a{t)q]. 
we may connect a temperature to this, 
aiTi 1 



T{t) 



At)' 



m 



ksTity 



(31) 



(32) 



In the approximation of instantaneous decoupling at tem- 
perature Ti the distribution function reads 



d^g E{q) 
(27rfi,)3 e^i"-'^ + 1 

\2irficJ J e^-'^ + l 



(36) 



ksTk 



The pressure is not taken from (23) but from (24), 



d^c 



(27rft)3 e'^'Jc-a + 1 2,E{q) 
{keTf r d^k fc^ 



(37) 



(27r^c)= 



+ 1 E{kBTk/c) 



This Ansatz has been adopted because it satisfles the 
energy conservation 



u + 3-{u + p) = 0. 



(38) 



Now our point is that this is consistent with thermody- 
namics when taking s = na from (28) and /U from (12), 



u+p-Ts , , q^c^ \ ^ 
n = \^(^^ + 3^/-^"' 



(39) 



f{q,t) 



1 



exp[/3i^(gTi/r)-ai] + l' 



(33) 



where fi\ = l/fc^Ti and ai = a(Ti) are time inde- 
pendent. For neutrinos the decoupling took place when 

they were relativistic, fcsTi ^ mc? . In that case we 
l3iE{qTi/T) f» (3qc, so f simplifies to 



/ 



1 



oPpc-a _j_ 2^ ' 



(34) 



where a = ai is a constant. Thus even though neutrinos 
have a mass, their density will be quasi-relativistic, at 
least as long as their distribution is uniform. (It has 
been argued that neutrinos are presently condensed on 
matter concentrations such as galaxy clusters [6].) 

It has long been supposed that in this non-equilibrium 
situation thermodynamics would still apply. 



where the averaging is over expression (33) or (34) for /. 
For ksT <gC mc^, where p/u « it yields /x « mc?, which 
is consistent with (11) because it cancels the large term 
u K, nmc?, thus repairing the relation (1). 

In conclusion, thermodynamics does explain the non- 
equilibrium situation of decoupled quantum matter in 
the expanding Universe, but the correspondence ji = 
aksT between the chemical potential /U and a, the time- 
independent parameter of the Fermi-Dirac distribution 
(34), is broken and replaced by the general relation (39). 
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